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INTRODUCTION 


Dynamical  stress  behaviors  and  shock  transients  in  mechanics  are  impor 
tant  subjects  to  be  studied.  The  use  of  the  finite  element  method  based  on 
developed  algorithms  from  the  variational  principle  can  give  direct  numerical 
solutions  for  partial  derivatives  of  the  functions  to  these  problems. 

Many  researchers  have  found  it  difficult  to  apply  the  finite  element 
method  to  hyperbolic  type  partial  differential  equations  (PDE).  Galerkin 
method  and  the  like  have  been  used  instead.  Previously,  the  author  has 
attempted  to  solve  these  hyperbolic  type  PDE's  by  employing  variational 
principles.  The  present  report  shows  that  under  certain  conditions,  the 
stationary  values  are  strong  extremals,  not  saddle  points.  This  is  equivalent 
to  requiring  that  the  second  variations  of  the  functional  be  positive  semi 
definite  with  discontinuities  in  the  partials. 

THE  FUNCTIONAL  BILINEAR  FORM 

A  functional  in  bilinear  form  (refs  1,2)  is  assumed  for  the  derivation  of 
the  Euler- Lagrange  equations  of  the  original  system  and  its  symmetrical 
adjoint  system  so  that  the  first  variations  of  the  functional  can  be  used  for 
numerical  computation  by  the  finite  element  method  (ref  3).  Let  us  consider 
the  functional 


lc.  N.  Shen,  "Method  of  Solution  for  Variational  Principle  Using  Bicubic 
Hermite  Polynomial,"  Transactions  of  the  27th  Conference  of  Army 
Mathematicians,  ARO  Report  82-1,  p.  247. 

2C.  N.  Shen,  "Variational  Principle  for  Gun  Dynamics  With  Adjoint  Variable 
Formulation,"  Proceedings  of  the  Third  US  Army  Symposium  on  Gun  Dynamics, 
Volume  II,  1982,  pp.  IV108-120. 

3A.  R.  Mitchell  and  R.  Wait,  The  Finite  Element  Method  in  Partial  Differential 
Equations ,  John  Wiley,  1977,  p.  143. 
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(1) 


J  ~  I  /  F[  t  ,x,y  ,yx,yj- ,yxx,y£t  ,y  ,yx,yt  ,yxx>yttl dxd^ 

to  xo 

where  y(x,t)  is  the  original  variable,  y(x,t)  is  the  adjoint  variable,  and  the 
subscripts  x  and  t  indicate  the  partial  derivatives. 

The  variations  of  the  variables  are  in  terras  of  a  small  parameter  y 
as 

y(x,t,y)  =  y(x,t)  +  <$y(x,t)  ,  <$y(x,t)  =  yn(x,t)  (2a) 

yx(x,t,  y)  =  yx(x,t)  +  «yx(x,t)  ,  «yx(x,t)  -  (2b) 

and 

yt(x,t, y)  =  yt(x,t)  +  6yt(x,t)  ,  6yt(x,t)  =  ynt(x,t)  (2c) 

Similar  expressions  can  be  obtained  for  higher  derivatives  of  y. 

The  functions  y(x,t,y)  and  y(x,t)  are  close  or  neighboring  in  the  sense 
of  closeness  (ref  4)  of  order  one  if  the  following  holds  for  some  e  >  0 

|yn(x,t)|  <  e  ,  lyr^x.t)!  <  e  >  and  |>JIlt(x*t)|  <  e  (3a) 

It  is  noted  from  above  that  the  partial  derivatives  of  the  variations  are 
small  for  a  weak  extremum. 

On  the  other  hand,  for  a  strong  extremum  we  require  only 

|yn(x,t)|  <  e  (3b) 

The  adjoint  variable  y(x,t,y)  and  its  partials  can  be  obtained  by  adding  a  bar 
on  top  of  the  original  variable  y  and  its  partials  as  given  by  Eq.  (2). 

THE  FIRST  VARIATION  AND  THE  EULER-LAGRANGE  EQUATIONS 

The  functional  J(y)  in  Eq.  (1)  can  be  expanded  in  a  Taylor  series  (ref  5) 
of  y  as 

J(y)  =  J(y=0)  +  6j  +  S2J  (4a) 

^L.  E.  Elsgolc,  Calculus  of  Variations,  Addison-Wesley ,  1962,  p.  15. 

^H.  Ruud,  The  Hamilton  Jacob  Theory  of  the  Calculus  of  Variations,  Robert  E. 
Krieger,  1973,  pp.  5-8. 
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where 


and 


<SJ 

SJ  =  p(--)p=o 

OU 


(4b) 


■2  32J 


62j  ■  r 


(4c) 


Equation  (1)  involves  two  independent  variables  x  and  t,  two  dependent 
variables  y  and  y,  and  each  variable  has  two  second  order  derivatives  with 
respect  to  x  and  t.  The  first  variation  can  be  deduced  from  Gelfand  and 
Formin  (ref  6)  and  separated  into  two  parts  as 

6J  =  6j(6y)  +  6j(<$y)  (5a) 

where 

tb  Xb  3  3  32  32 

Wy)  -  J  I  (Fy  -  --  Fyx  -  --  Fyx  +  Fytt  +  2  Fytt 

tQ  xQ  3t  9x  St'" 


3t  3x 


and 


+  ~~2  Fyxx) 6y(x,t)dxdt 

oX 


tb  xb  -  3  -  3  -  32  -  32  - 

6j(6y)  =  j  j  (Fy  -  --  Fyt  -  --  Fyx  +  Fytt  +  2  —  Fytx 
tQ  xQ  3t  3x  3t  3t  3x 


(5b) 


+  ~2  Fyxx)6y(x,t)dxdt 
ox 


(5c) 


For  the  wave  equation  let  us  consider  the  following  F  in  Eq.  (1) 


1  -  7  1  - 

f  =  -  y(ytt~a  yxx)  +  yQ  +  -  y(ytt_a  yxx)  +  yQ 


(6) 


By  evaluating  the  terms  in  Eq.  (5b)  for  the  above  F#  we  have 

tb  fxb  - 

6j(  6y)  =  J  J  <f>  <5y  dxdt 
to  xo 


(7a) 


^1.  M.  Gelfand  and  S.  V.  Formin,  Calculus  of  Variations,  Prentice-Hall,  1963, 
pp.  22,  34,  42,  and  161. 
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where 


(7b) 


* ■  Ytt  ■  a2yxx  +  Q 


Similarly  by  Eq.  (5c)  one  obtains 


(7c) 


where 


♦  ■  y«  -  a2yxx  +  Q 


(7  d) 


The  first  variations  of  6y  and  <5y  are  arbitrary  within  certain  limitations  to 
be  discussed  later.  Then  the  fundamental  lemma  (ref  7)  of  calculus  of 
variation  gives 


♦  -  ytt  ”  a2yxx  +  Q  =  0 


(8a) 


and 


^  =  vtt  ~  a2yxx  +  Q  =  o 


which  are  the  adjoint  and  the  original  systems  of  the  wave  equation.  Note 
that  the  above  two  equations  are  similar  and  interchangeable  by  adding  or 
dropping  the  bar  on  top  of  y. 

THE  SECOND  VARIATIONS 

In  order  that  the  functional  J  in  Eq.  (1)  be  an  extremum  (ref  8), 
the  second  variation  <5^J  must  be  either  positive  semi-definite  for  a  minimum, 
(or  negative  semi-definite  for  a  maximum),  i.e.. 


S2J  >  0  (or  62J  <  0) 


(9) 


The  above  is  the  second  necessary  condition  for  a  minimum  (or  a  maximum). 


^R.  Weinstock,  Calculus  of  Variations,  McGraw-Hill,  1952,  p.  16. 

®Hans  Sagan,  Introduction  to  the  Calculus  of  Variations,  McGraw-Hill,  1969, 
p.  38. 
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From  Eqs.  (4b),  (7a),  and  (7c),  one  obtains 


Sj  fcb  xb  - 

—  =  J  j  ( <}>  6n  +  4>  <5n)dxdt  (10) 

9  y  tQ  xQ 

where 

<t>  =  Q  +  yntt  -  a2yTlxX  (Ha) 

and  -  -  -  - 

<j>  =  Q  +  yntt  -  a2yn*x  (lib) 

92J  tb  xb  9x  9y  9x  9yt  9x  ^x 

9y2  t0  xQ  ^ 9y  9y  9yt  9y  9yx  9y 

9x  8ytt  9x  9ytx  9x  9yxx  - 

+ - +  2 - + - ]  dn 

9ytt  9y  9yxx  9y  3yxx  3p 

+  (similar  term  with  bars)dh}dxdt  (12) 


—  j  j  ^(  ^tt- a^r'xx)  0  4"  (  H^-a  2nxx)  nldxdt  (13) 

to  x0 


Integrated  by  parts  we  have 


2  92J  ,tb  xb 


1 


tQ  xr 


y2[-nt(x,t)rit(x,t)  +  a2nx(x,t)  nx(x,t)]dxdt 
+  B.C.  +  I.C. 


(14) 


It  can  be  proved  that  for  physical  problems  such  as  a  fixed  or  a  free  end  on 
either  side,  the  boundary  conditions  are  zero.  The  initial  conditions  are 
zero  as  illustrated  in  the  next  section. 


THE  ADJOINT  SYSTEM 

The  adjoint  system  may  be  taken  as  the  image  reflection  in  the  time 
domain  of  the  original  system,  as  shown  in  Figure  1. 


y(x,t)  =  y(x,tb+t0-t) 

(15a) 

yx(x>t)  =  yx(x>tb+tQ-t) 

(15b) 

yt(x,t)  =  -yt(x,tb+t0-t) 

(15c) 
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In  addition,  Eq.  (15)  yields  the  known  initial  conditions  as 

y(x,tb)  =  y(x,tQ)  (16a) 

yt(x,tb)  =  -yt(x,t0)  (16b) 

Since  the  adjoint  system  is  a  reflected  mirror  in  time,  the  far  end  initial 
conditions  for  the  adjoint  system  are  known. 

We  may  now  derive  the  variations  of  Eqs.  (15)  and  (16),  which  give 


Sy(x,t)  =  5y(x,tb+t0-t) 

(17a) 

fyxU.t)  =  6yx(x>tb+to"t) 

(17b) 

<5yt(x,t)  =  -<5yt(x,tb+t0-t) 

(17c) 

<5y(x,tb)  -  6y(x,t0)  =  0  for  all  x 

(18a) 

6yx(x,tb)  =  -6yt(x,t0)  =  0  for  all  x 

(18b) 

By  substituting  Eq.  (15)  into  Eq.  (14)  and  using  Eq.  (2),  we  have 


where 


62J  = 


*b 

J  P(x,t)dxdt 


P(x,t)  -  <5yt(x,t)6yt(x,tb+t0-t)  +  a26yx(x,t)  6yx(x,tb+t0-t) 
=  p2[nt(x,t)nt(x,tb+t0-t)  +  a2nx(x,t) ^(x.th+tQ-t)] 


(19a) 


(19b) 


SENSITIVITY  RELATIONSHIP 

In  order  to  show  that  the  second  variation  of  the  functional  J  is 
positive  semi- def inite,  one  needs  to  obtain  the  variations  of  the  function  and 
its  partials  together  with  those  of  the  adjoint  function  and  its  partials  as 
indicated  in  Eq.  (19).  We  can  get  these  variations  through  the  study  of  the 
sensitivity  coefficients  (ref  9)  and  their  relationship  to  the  parameters 


^Rajko  Tomovic,  Sensitivity  Analysis  of  Dynamic  Systems,  McGraw-Hill,  1963. 
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given  in  Eq.  (2).  Let  the  forcing  function  in  Eq.  (7d)  be 

Q(x,t)  =  qf(x,t)  (20) 

It  is  assumed  that  the  forcing  function  parameter  q  is  subject  to  a  small 
constant  perturbation  <5q  as 

q  =  q0  +  6q  (21) 

Then  the  variation  of  the  function  y  is 


where 


«y(x,t) 


9y(x,t) 

3q 


6q  =  v(x,t)  <$q 


3y 

v(x,t)  =  - 

3q 


(22a) 


(22b) 


The  quantity  v  is  the  sensitivity  coefficient  for  the  variation  6y(x,t)  due  to 
a  small  constant  perturbation  Sq. 

The  original  PDE  in  Eq.  (7d)  can  be  written  as 
4>  =  Ly  +  Q 

=  ytt  "  a2yxx  +  qf(x»t)  ■  0  (23) 

Due  to  the  perturbation  of  q,  the  change  of  <j>  obeys  the  following  relation¬ 
ship: 

9<t>  3Ytt  3  4>  9yxx 

- - 4- - -  +  f  (x,  t)  =  0  (24) 

3ytt  9q  9yxx  9q 

It  is  also  noted  from  Eq*  (23)  that 


3  <j) 

3(j> 

-  =  1  and 

— — 

9ytt 

9yxx 

Using  the  definition  in  Eq#  (22b),  the  partials  can  be  interchanged  as 


9ytt 

9q 


32  3y 

T“  (“)  =  Vtt 

9q 


3t2 


9yXX  3 2  3y 

and  -  — -5  ( — )  =  vtt 

8q  3x2  3q  C 


(25) 


(26) 
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Substituting  Eqs.  (25)  and  (26)  into  Eq.  (24),  we  have 

vtt  -  a2vxx  -I-  f(x,t)  =  0  (27) 

If  we  compare  the  definitions  of  variation  in  Eq,  (2a)  with  the  definition  of 
sensitivity  relationship  in  Eq.  (22a),  we  have 

Sy(x,t)  =  un(x,t)  -  (<$q)v(x,t)  (28) 


which  gives 

n(x,t)  =  v(x , t )  (29a) 

and 

Sq  =  v  (29b) 

Thus  Eq.  (27)  becomes 

ntt  -  a2nxx  -4-  f(x,t)  =  0  (30) 

which  gives  the  PDE  of  the  variations  of  the  original  system. 

If  we  compare  Eq.  (30)  with  Eq.  (23),  we  see  that  the  variation  n(x,t)  * 
y-16y(x,t)  in  Eq.  (30)  takes  the  place  of  the  function  y  in  Eq.  (23)  with  q  ■ 
1.  Therefore,  the  PDE  for  the  variations  is  unchanged  except  by  a  scale 
factor.  Thus  the  solution  of  the  variation  <$y(x,t)  has  the  same  form  as  that 
of  the  original  function  y. 

Similarly  for  the  adjoint  system,  one  can  obtain 


Sy(x,t)  =  un(x,t)  =  (<5q)v(x,t) 

(31) 

n(x,t)  =  v(x,t) 

(32a) 

6q  =  u 

(32b) 

Htt  "  a2nxx  +  =  0 

(33) 

which  is  the  PDE  of  the  variations  of  the  adjoint  system. 
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THE  WEAK  EXTREMUM 


Let  us  assume  that  the  variations  are  separable  in  space  x  and  time  t, 

i.e. , 


n(x,t)  =  f(x)g(t)  n(x,tb+tQ-t)  =  f(x)g(tb+tQ-t)  (34a) 

nt(x,t)  =  f(x)gt(t)  nt(x,tb+t0-t)  =  f(x)gt(tb+t0-t)  (34b) 

nx(x,t)  =  fx(x)g(t)  nx(x,tb+t0-t)  =  fx(x)g(tb+t0-t)  (34c) 

Substituting  Eq.  (34)  into  Eq.  (19b)  gives 

P(x,t)  =  f2(x)gt(t)gt(tb+t0-t)  +  a2fx2(x)g(t)g(tb+t0-t)  (35) 

with  the  initial  conditions  from  Eqs.  (18a)  and  (18b) 

n(tb)  =  n(t0)  =  0  =>  g(t0)  =  0  (36a) 

nt(tb)  -  -nt(t0)  =  o  =>  gt(t0)  =  0  (36b) 

Let  us  construct  the  function  g(t)  with  tQ  =  0,  and  w  >  0,  as 

g(t)  =  1  -  cos  wt  ,  g(tQ)  =  0  (37a) 

and 

gt(t)  =  0)  Sin  Wt  ,  gt(tQ)  =  0  (37b) 

which  satisfy  Eq.  (23).  Then  we  have 

gt(t)  >0  if  0  <  ait  <  wtj-j  <  tt  (38a) 

gt(tb-t)  >  0  if  0  <  u<tb-t)  <  wtb  <  *  (38b) 

g(t)  >0  if  0  <  (ot  <  wtb  <  2 tt  (39a) 

and 

g(tb-t)  >  0  if  0  <  a<tb-t)  <  a>tb  <  2  tt  (39b) 

Under  these  conditions  the  second  variation  is  positive  semi-definite,  i.e., 

P(x,t)  >  0  ,  if  <tb-t0)  <  TT  (40a) 

or  tb  -  tQ  <  tt/o)  (40b) 


The  above  g^(t)  is  continuous  and  requires  close  neighborhood  of  the 
function  y^.  This  is  a  weak  extremum.  However,  actual  solutions  of  the  wave 
equation  show  there  are  jumps  in  the  functions  yt  and  yx.  This  can  only  be 
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proved  by  a  strong  extremum  later  in  the  report.  Let  us  assume  in  Eq.  (34a) 
that 

f(x)  =  (xjj-x)  2(x-x0) 2  >  0  for  xQ  <  x  <  xb  (41) 

which  satisfies  all  possible  boundary  conditions  for  the  variations  of  the 
wave  equation.  It  is  also  noted  that  for  the  chosen  g(t)  in  Eq.  (37),  the 
variation  6y(x,t)  in  Eq.  (7a)  and  its  partials  obtained  from  Eqs.  (2a)  and 
(34a)  are  always  arbitrarily  positive  for  tb  <  n/u>.  This  also  satisfies  the 
fundamental  lemma  of  the  calculus  of  variations  which  gives  the  Euler-Lagrange 
Eq.  (8a).  A  similar  approach  follows  that  Sy(x,t)  >0  in  Eq.  (7c)  for  the 
derivation  of  Eq.  (8b). 

PERIOD  OF  THE  WAVE  EQUATION 

It  is  not  true  that,  "For  a  vibrating  string  there  is  no  time  interval 
short  enough  to  guarantee  that  y(x,t)  actually  minimizes  the  action 
functional,"  (ref  6).  This  conjecture  was  based  on  the  series  solutions  of 
the  wave  equation,  as  a  system  of  infinitely  many  coupled  oscillators,  with 
infinitely  many  natural  frequencies,  (i.e.,  m  •>  “) .  However,  if  the  problem 
is  approached  from  a  different  angle,  the  solution  of  the  wave  equation  with 
no  forcing  terms  is  given  as  the  sum  of  two  wave  fronts  (ref  10): 

y(x,t)  =  X(x-at)  +  p(x+at)  ,  for  Q  =  0 
y(x,0)  =  X(x)  +  p(x)  =  <|>(x) 

=  -aXx(x)  +  apx(x)  =  3(x) 

x-0 

^1.  M.  Gelfand  and  S.  V.  Formin,  Calculus  of  Variations,  Prentice-Hall,  1963, 
pp.  22,  34,  42,  and  161. 

R.  Wylie,  Advanced  Engineering  Mathematics,  McGraw-Hill,  1951,  pp.  211— 

212. 
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Then 


(42) 


1  1  j&at 

y(x,t)  =  -  [  4>(x-at)  +  4>(x+at)]  +  —  J  6(s)ds 

2  2a  x-at 

Let  us  take  an  example  of  a  prismatic  bar  (ref  11)  with  the  left  end  fixed  and 

the  right  free.  If  the  free  bar  is  initially  stretched  by  a  constant  force 

and  suddenly  released,  the  stress  wave  is  travelling  from  right  to  left  at  a 

period  T  *  4£/a,  where  &  =  x^  -  xQ  is  the  length  of  the  bar.  From  Figure  2 

one  can  write  in  terms  of  unit  step  functions  u  for  the  partial  derivatives  yx 

and  yt>  for  the  wave  fronts  at  some  intermediate  location  K  and  time  5. 

yx(x,t=5)  =  [u(x-O)  -  u(x- K)]  (43a) 

T 


yt(x,tss5)  =  -a(u(t-0  -  u [  t  -  -  £)]}  (43b) 

where  the  period  T  =  4£/a  (43c) 

Z/l  =  |  (T/4)  -  c|  / (T/4)  (43d) 


The  corresponding  values  of  5  and  5  are 

at  5  =  £,  £/2,  0,  -2-/2,  and  & 


(44) 


Then  T/8 ,  T/4,  3T/8,  and  T/2  respectively. 

It  is  noted  that  the  "physical  frequency"  of  the  system  is 

(op  -  2tt/t  =  (Tra)/(2£) 

This  frequency  is  a  constant  never  approaching  infinite  for  finite  parameters 
a  and  £. 


THE  STRONG  EXTREMUM 

For  the  variations  of  Eq.  (43),  shown  in  Figures  2  and  3,  we  let 

6yx(x>0  =  yx(x>  -  yx(x»0 

=  y[u(x-0)  -  u(x-£)]  (45) 


H-L.  S.  Jacobsen  and  R. 
473  (Figures  10-18). 


S.  Ayre,  Engineering  Vibrations,  McGraw-Hill,  1958,  p. 
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for  0  <  5  <  i9 


<5yx(x,0  “  0  at  the  right  end 
<5y(x,C)  =0  at  the  left  end 


Thus  Eq.  (45)  satisfies  the  required  boundary  conditions.  Similarly, 

<$yt(x,  O  "  -ay{u(t-C)  -  u[  t-T/2)  ]  } 


(46) 


for  C  >  0, 


6yt(x,C)  =  0  and  6y(x,C)  ~  0  initially 


It  is  noted  that  both  6yx  and  <Syt  are  constant  within  the  intervals.  Thus  Eq. 
(46)  satisfies  the  required  initial  conditions.  Equations  (45)  and  (46) 
fulfill  the  “general  conditions”  for  the  variation  and  are  subject  to 
discontinuity  of  the  partial  derivatives,  where  the  step  functions  occur. 

This  leads  to  a  strong  extremum. 

We  use  the  increment  in  time  for  computation  by  the  finite  element  method 


as 


tb  ~  tQ  =  a  »  nT/4 


(47) 


where  n  *  1  or  2. 

Under  this  condition  if  t  varies  from  t0  to  (t0+a),  then  a  -  t  varies  from 
(to+tf)  to  tQ.  We  further  divide  the  increment  into  four  equal  intervals,  each 
interval  having  T/16  in  time. 

Since  <5yx  and  6yt  are  constants  or  zeroes,  the  second  variation  as  given 
by  Eq.  (19)  can  be  changed  into  finite  summation  instead  of  integration  over 
the  time  domain  as, 


xb  N-l  xb 

(62j)k  =  ( At/2) /  Pk(x,t)dx  +  (At)  l  j  Pk+i(x,t)dx 
xQ  1=1  xQ 


,xb 

+  (At/2)  /  Pk+N(x,t)dt 
xo 


(48) 
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where  k  is  the  index  of  advance  and  N  is  the  number  of  intervals.  For  the 
interval  of  T/16  and  with  the  aid  of  Eq.  (47),  the  quantity  N  becomes  N  =  4  if 
n  =  1  and  N  =  8  if  n  =  2. 

From  Eq.  (19b)  we  have 

Pk+i(x,t)  =  6yt[x,(t0+kT/16+iT/16)]6yt[x,(t0+kT/16+o-iT/16)] 

+  a26yx[x,(t0+kT/16+iT/16)] 6yx[x,(to+kT/16+0-iT/16)] 

A  (49) 

for  i  =  0,1,2,3,...N 
k  =  0,1,2,3, . . . 

If  we  integrate  the  above  products  in  the  spatial  domain  for  n  =  1,  we 
obtain  the  the  following  as  shown  in  Figures  4(a)  and  4(b). 


,xb 

J  Pk(x,t)dx  =  0  for  all  k 
xo 
,xb 

J  Pkq.4(x,t)dx  =  0  for  all  k 
xo 
xb 

/  Pk+i(x,t)dx  =  a2y2  A/2  >  0  for  all  k 


xb 

J  Pk+3(x,t)dx  =  a2y2  V2  >  0  for  all  k 


(50a) 

(50b) 

(50c) 


(50d) 


and  xb 

/  Pk+2(x>t)dx  =  a2U2£  >  0  for  a11  k  (50e) 

xo 

Substituting  Eq.  (49)  into  Eq.  (48)  gives 

(Sj2)n=i  =  ( AT)  [0  +  a2&/2  +  a2Jl/2]  y2(  At)  >  0  (51) 

For  n  =  1,  i.e.,  the  increment  a  being  (T/4)  the  quarter  period,  the  second 
variation  is  positive  semi- def inite.  Thus,  the  variational  method  is  valid 
and  the  functional  yields  a  minimum. 
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The  integration  of  Eq.  (49)  in  the  spatial  domain  for  n  =  2  (i.e.,  o  = 


T/2)  gives  the  following  for  all  k  as  shown  in  Figures  5(a)  and  5(b). 


xb 

/  Pk(x,t)dx  =  -a2w2& 
xo 
xb 

/  Pk+i(x,t)dx  =  -a2P2A/2 
xo 

xb 

J  Pk+2(x»t)dx  =  0 


xb 

/  Pk+3(x,t)dx  =  a2&P2/2 
xo 


(52a) 

(52b) 

(52c) 

(52d) 


and 


Pk+4(x,t)dx  =  a2^y2/2 


for  all  k 


(52e) 


Substituting  Eq*  (52)  into  Eq .  (47)  gives 

( 62J)n=2  =  0  for  all  k  (53) 

For  n  =  2,  i.e.,  the  increment  o  being  (T/2)  the  half  period,  the  second 
variation  is  zero.  Thus  the  variational  method  is  on  the  margin  to  yield  a 
minimum  for  the  functional.  It  can  be  proved  that  for  o  >  T/2  the  second 
variation  becomes  negative.  Thus  the  variational  method  is  valid  only  if  the 
increment  o  used  for  the  finite  element  method  is 


o  <  T/2 


(54) 


where  T  is  the  period  corresponding  to  the  "physical  frequency"  of  the  system. 


CONCLUSIONS 

The  functional  in  bilinear  form  is  established  and  the  Euler-Lagrange 
equations  are  derived  to  obtain  the  wave  equation  and  its  symmetrical  adjoint 
equation.  The  second  variation  of  the  functional  is  given  in  terras  of  the 
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variations  of  the  partials  in  space  and  time,  both  in  the  original  variable 
and  the  adjoint  variable.  The  adjoint  system  is  defined  as  the  image 
reflection  in  the  time  domain  of  the  original  system.  The  physical  frequency 
of  the  system  is  determined  (not  series  expansion  frequency).  There  is  an 
upper  bound  for  the  increment  in  time  that  can  be  used  for  computation 
purposes.  The  strong  extremum  was  illustrated  by  an  example  of  the  wave 
equation  where  the  first  partials  of  the  variables  and  their  variations  were 
subject  to  discontinuities. 
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Figure  1.  Image  Reflection  of  the  Adjoint  System. 
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Figure  2.  Spatial -Temporal  Response  of  a  Prismatic  Bar  With 
Sudden  Release  of  Load  at  One  End. 
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Figure  3.  Spatial  and  Time  Derivatives  for  Figure  2. 
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Figure  A.  Integral  of  Products  of  Spatial  and  Time  Derivatives 

Using  Quarter  Period  Increment  for  Finite  Element  Method. 
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